Energy of gravitational radiation in plane-symmetric space-times 
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Gravitational radiation in plane-symmetric space-times can be encoded in a complex potential, 
satisfying a non-linear wave equation. An effective energy tensor for the radiation is given, taking 
a scalar-field form in terms of the potential, entering the field equations in the same way as the 
matter energy tensor. It reduces to the Isaacson energy tensor in the linearized, high-frequency 
approximation. An energy conservation equation is derived for a quasi-local energy, essentially the 
Hawking energy. A transverse pressure exerted by interacting low-frequency gravitational radiation 
is predicted. 
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I. INTRODUCTION 



Gravitational radiation, as predicted by Einstein gravity, is indirectly observed in such examples as the Hulsc- Taylor 
pulsar, and widely expected to be directly observed in the coming years, offering a new window to understand various 
astrophysical processes, such as binary inspiral and merger of black holes or neutron stars. However, the textbook 
theory of gravitational radiation mostly concerns weak radiation, either in the linearized approximation or at infinity 
in an asymptotically flat space-time 1]. Comparatively little is known about strong- field radiation. One exception 
is plane gravitational radiation, where exact solutions describe radiation propagating in one direction. The simplest 
qh' scenario to study interaction effects is the head-on collision of two such beams, as pioneered by Szekeres 0, 0] and 
1 , reviewed by Griffiths pj. More generally, one may study plane symmetric space-times, which in vacuum generally 
^jj^ consist of gravitational radiation propagating in opposite directions and interacting [5] . 

, Much is known about such space-times, including that the interaction is non-linear, that the key dynamical equations 
can be cast as a complex Ernst equation [6i], and that the cross-focusing of the radiation produces a caustic which 
is generically a curvature singularity, though there are non-generic exceptions 0, Q ■ This article introduces an 
effective energy tensor Q for the gravitational radiation, taking a scalar-field form in terms of a complex potential 
<I>. Then O enters the field equations in the same way as the matter energy tensor, in particular entering an energy 
QQ , conservation law. The Ernst equation is manifestly a wave equation for $, generally with a non-linear source, which 
^SJ ' vanishes for coUinear polarization. 

jy-^ , The method involves a conserved time vector /c°, a conserved energy-momentum density j", a corresponding energy 
' E and a first law for E involving energy-supply and work terms. Surface gravity k is also defined and takes a quasi- 
00 ' Newtonian form. This is intended to complete the same pro gram me of identifying physical quantities and equations 
which has previously been performed in spherical symmetry [lOl . |Tl| , cylindrical symmetry and a quasi-spherical 
' approximation [l^, O, [l^, [3| • These references will be assumed for comparison throughout the text without repeated 
citation, though the treatment here is self-contained. 



II. METRIC VARIABLES AND FIELD EQUATIONS 



Cartesian coordinates (z, y) on the planes of symmetry will be used, to allow easy comparisons with standard 
coordinates (z, (p) in cylindrical symmetry and Lp) in spherical symmetry and the quasi-spherical approximation, 
ft is convenient to use null coordinates in the normal space, as they are adapted to gravitational radiation. Then 
the metric can be written locally as 

ds^ = -^e^"' dx+ dx~ + A (e^* sec 2xdy'^ + 2 tan 2xdydz -)- e"^"^ sec 2xdz'^) (1) 

where (A, 0, x,7) are functions of (a;+,a;^). Here A is the specific area, meaning that it is the area of a square 
coordinate patch (0, 1) x (0, 1) in the {y, z) plane. It is invariant up to constant linear transformations of y and 
z, under which it scales by a constant factor. The remaining freedom in (j/, z) is by rotations, under which A is 
invariant. The functions (0, x) encode the gravitational radiation, as will be seen below. They are invariant up to 
the above-mentioned transformations of {y,z), which will be treated as fixed henceforth. The remaining function 7 
is invariant up to functional rescalings x^ 1— > x^ {x^ ) , under which it transforms by additive functions of x^ and x~ . 
The variables have been chosen so that the induced metric on the planes of symmetry takes a similar form to that 
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used in the quasi-spherical approximation, with {dz,dy) replaced by [M^sm'ddLp), and takes a similar form to that 
used in cylindrical symmetry. The Szekeres variables (P, A/, Q, W) are related by 

P = - log A, M = -27, Q = -20, sinh W = tan 2x (2) 

or cosh W — sec 2%. 

The six independent components of the Einstein equation may be found directly, or by comparison with the Szekeres 
form, as 

2Ad±d±A ~ {d±Af - AAd±Ad±^ + AA^ sec^ 2x {{d±(t>f + {d±xf) = -1&ttA^T±± (3) 

d+d-A^ %ttAT+^ (4) 

2Ad+d-4> + d+Ad-(t> + a_Ac»+?!) + AAt&u2x{d+(t>d-x + = AnAe^'^ cos^ 2x(T^^ - P;f) (5) 

2Ad+d-X + 9+^a_x + d-Ad+x + 4Atan2x(a+x9-X - = AttAc^'^ cos^ 2x{e^'^Ty + e-^-^T^^) (6) 

4v42a^a_-^ _ + AA^ sec^ 2x(9+(/'9_0 + d+xd-x) = -Stt^I^ (2T+_ + e^'^iT^ + T^^)) (7) 

where d± = d/dx^, T denotes the energy tensor of the matter with T±± — T{d±, d±), = P(9+, and the units 

are such that Newton's gravitational constant is unity. The equations ([3]) can be regarded as constraint equations on 
initial null hypersurfaces S± of constant x^, as they are preserved in the i9ip directions due to the Bianchi identities 
or energy- momentum conservation. The other equations (H])-® are then the evolution equations. 

In vacuum, T = 0, it is well known that these equations describe the propagation and interaction of gravitational 
radiation in the opposite d± directions, and that the radiation may be encoded in (0, x)- The solution to ([4]) is trivial 
and can be used to fix the rescaling freedom in x^. One may give initial data for {(j),x) on T,±, corresponding to 
initial radiation profiles, with ([3]) determining 7 on Then the main task is to solve (I5|)-(l6|) simultaneously for 
(0, x), after which the full solution follows from ([7]) by quadrature for 7. The main equations ©-(IH]) can be written 
as a complex Ernst equation, corresponding physically to a non-linear wave equation, as will be verified below. 



III. EFFECTIVE ENERGY TENSOR FOR GRAVITATIONAL RADIATION 



The next aim is to find an effective energy tensor Q for the gravitational radiation, analogous to those found in 
cylindrical symmetry and the quasi-spherical approximation, and consistent with the Isaacson effective energy tensor 
in the high-frequency linearized approximation In all cases, the components of the energy tensor are quadratic in 
first derivatives of the metric, in this case the d± derivatives of {4>,x)t and such terms can be seen in the last term 
in parentheses on the left-hand side of each of ([3]), ([S])-©. The idea is to identify these terms as components of the 
desired Q, corresponding to the components of T on the right-hand sides. The result is that one may introduce a 
complex potential 



<^ = (p + iX (8) 



and define the effective energy tensor as 



2V(,$Vfc)$ - g,fa5=''Ve$Vrf$ 

"ab = T = ■ (y) 

87rcosh2($- $) ^ ' 

where g is the space-time metric and V its covariant derivative operator. It is manifestly a tensor, taking a scalar-field 
form in terms of (f>, with the same form, including the same denominator, as in the quasi-spherical approximation. 
Apart from this denominator, it is the energy tensor of a massless complex scalar field $. Explicitly in terms of ((/>, x), 

2Va0Vb0 + 2VaX^bX - 9abg^'^{yc(bVd^ + V.xVdX) ^^ 

Bah = 5 T7i ■ (10) 

Sttcos^ 2x 

If X = 0, it reduces to the energy tensor of a massless scalar field (f>, as in cylindrical symmetry, where the corresponding 
(j) reduces to the Newtonian gravitational potential in the Newtonian limit. Here there are generally two polarizations 
of the radiation, as is familiar from the linearized approximation. Inspection of the metric (H]) for small $ identifies (j) 
as encoding the "plus" polarization and x as encoding the "cross" polarization. These properties justify the numerical 
factors chosen in the definitions of (0, x) and partly motivated the chosen symbols. 
The non-trivial components of follow explicitly as 

4^e±±=sec2 2x((9±(/.)' + (9±x)') (11) 

e+_ = o (12) 

47r_Le = e-^'' sec^ 2xid+(j)d-4> + d+xd-x)^9 (13) 
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where _L denotes projection onto the planes of symmetry and the transverse metric is given in (y, z) coordinates by 

It is then straightforward to verify that adding to T on the right-hand sides of the Einstein equations ©-([T]) 
cancels the quadratic terms in (0, x) on the left-hand sides. In abstract terms, the Einstein equation G = SirT may 
be rewritten as C = 87r(T -I- O) in terms of a truncated Einstein tensor C, whose components have a simpler form to 
those of the Einstein tensor G. 

The physical interpretation of &±±/2 is the energy density of gravitational radiation propagating in the direction. 
Apart from the non- linear modification due to the sec^ 2% factor in (fTTj) . it is the energy density of a complex scalar 
field The numerical factor also corresponds to the energy density of electromagnetic radiation in Gaussian units, 

with (p corresponding to the electric potential and x vanishing. The vanishing of (|12p is familiar from cylindrical 

symmetry and the quasi-spherical approximation, and indicates that the gravitational radiation is workless. Note that 
this is generally not so for a similar effective energy tensor found in the context of black holes [13, [3 ^^'^ uniformly 
expanding flows [l^, [l^. The non-negativity of 8±± indicates that, as an energy tensor, Q satisfies the dominant 
energy condition, meaning physically that gravitational radiation carries positive energy. The other non-zero terms 
(|13p indicate that interacting gravitational radiation generally exerts transverse pressure and shear, proportional to 
the transverse metric. These terms vanish for radiation propagating in one direction only, where <I> is a function of 
a;"*" (or x~) only. They are commonly known as plane waves, but since this would appear to imply periodicity in some 
sense, this article uses the more general terminology of radiation. 



IV. CONSERVATION OF ENERGY 



To see how Q further qualifies as an effective energy tensor, one may proceed by analogy with spherical symmetry, 
cylindrical symmetry and the quasi-spherical approximation. Here the definitions and equations will be stated first 
in a manifestly invariant way, then verified in coordinates. First introduce the specific area radius 



r = \AV4^. (15) 

This is defined in order to compare with spherically symmetric space-times or the quasi-spherical approximation, so 
that one may easily treat astrophysical gravitational radiation as observed on or near Earth, since distant sources can 
be treated as points, producing roughly spherical wavefronts which can be treated as planes when observed. 
The Hodge operator * defines the Hodge dual *a of a normal one-form, up to sign, by 

g~^{*a,a) = 0, g^'^{*a,*a) ^ -g^^{a,a). (16) 

Then a preferred time vector is defined by 

k^g-^{^dr) (17) 
where the qualification "specific" is omitted here and henceforth. This vector is conserved: 

V-A; = 0. (18) 

The corresponding energy-momentum density is 

j = -g~\{T + e)-k). (19) 

Then j is also conserved: 

V • j = 0. (20) 

Here the standard physical interpretation is conservation of energy, and the role of O as an effective energy tensor is 
clear in that it appears additively with T in j. 

Put another way, both k and j are Noether currents, and the corresponding Noether charges are area volume 

F=|7rr3 (21) 

and energy defining the latter. Specifically: 

Ag{k) = *dV, Ag{j) = *dE. (22) 
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Integrating for E and requiring it to vanish for fiat space-time, 

E = -^rg-\dr,dr) (23) 

which has a similar form to the Misner-Sharp energy in spherical symmetry and the modified Thorne energy in 
cylindrical symmetry. In fact, if the planes of symmetry are toroidally compacted by periodic identifications in (y, z) 
at and 1, so that A is the area, then E coincides with the Hawking energy [2l| . 

Note that E > for trapped surfaces, E — for marginal surfaces and £^ < for untrapped surfaces. In particular, 
E vanishes for radiation propagating in one direction only. Thus it should not be interpreted as the energy of a wave 
in any sense. Taking the example of two colliding beams, where the surfaces in the interaction region are trapped 
if the null energy condition holds, one may interpret E as measuring energy due to cross-focusing of radiation. In 
particular, it diverges at the caustic formed by such cross-focusing. 

Introduce the work density 

u; = -trr/2 (24) 

and the energy flux 

tp = {T + Q) ■ g-\dr) + wdr (25) 

where the trace is in the normal space. Then conservation of energy (|20p can be written in the form of a first law: 

dE = Ail) + wdV (26) 

which has the same form as in spherical symmetry and the quasi-spherical approximation. Here the two terms can be 
interpreted as energy supply and work respectively, as in the first law of thermodynamics. Note again that Q appears 
additively with T in ip and (in a null sense) w, playing the role of an effective energy tensor. 
The corresponding definition of surface gravity is 

K ^ *d*dr/2 (27) 

where d is the exterior derivative of the normal space. Then the Einstein equations yield 

E , , 

K = — - 47rrw (28) 

which again has the same form as that in spherical symmetry and the quasi-spherical approximation. Apart from the 
matter term, this has the form of Newtonian gravitational acceleration. 
In dual-null coordinates ([T]), the corresponding expressions are 

— —O-j^dx^ + a^dx'^ where a — oij^dx'^ + a^dx~ (29) 

k = e-'^^d+rd- - d-rd+) (30) 

J ^ e-47[((T__ + e__)a+r - T+_a_r)5+ - ((T++ + e++)5_r - T+^d+r)dJ\ (31) 

E = e-'^^rd+rd-r (32) 

w ^ e"2''T+_ (33) 

V'i - -e-'^(T±± + Q±±)d^r (34) 

K = -e-'^^d+d-r. (35) 

Writing 4:{4:n)^/^E — e""^^ A^^^'^d+Ad-A and using the Einstein equations ©-(HI), a calculation yields 

d±E = Ae-^\d±rT+^ - d^r{T±± + e±±)). (36) 

Comparison with 

Ag{j) = Ae-^''[{{T++ + e++)d-r - T+^d+r)dx+ - ((r__ + e__)9+r - T+_9_r)d2;-] 

= [-d+Edx+ + d-Edx-] = 't'dE (37) 
verifies (|22p . Similarly, the calculation 

A{yj± + wd±r) = Ae-^\ - d^r{T±± + e±±) + 9±rT+„) = d±E (38) 
verifies (|26p . The easiest way to verify the conservation equations (fT5|) . (|20p is to use and exterior calculus: 

V • fc = A-^*d*{Ag{k)) = A-^^d^ifdV = (39) 

V • j = A-^*d*{Ag{j)) = A-^*d**dE = (40) 

since = ±1 and dd = 0. Finally, a calculation using the Einstein equation ^ verifies 



5 



V. GRAVITATIONAL WAVE EQUATION 

As is well known, the propagation equations (O-® for {(f>, x) can be written as a single complex Ernst equation, 
usually given in terms of an Ernst potential Z — e^* ot E = tanh <f> 4] . The corresponding form for $ is 

= 2tanh($- $)g"^(V$, V$) (41) 

where _LT = for simplicity. This has the same form as that in the quasi-spherical approximation. The calculation 
is straightforward: 

= (2a+a_$ + + (42) 

and 

2tanh($- $)5"i(V$,V$) = -Ap-^^ t£inh2ixid+(t) + id+x)id-<j) + id-x) 

= 4e-27 tan2x i{d+<pd.x + d-C^d+x) + i{d+xd-X - d+<pd-<P)) (43) 

then the result follows by comparing with (O-®. 

Note that (|4T|) is manifestly a wave equation for <I>, equating to a non-linear term in $. This source term is 
highly non-linear, being quadratic in V$ and also involving tanh(<I> — $). In the special case of coUinear polarization 
X = 0, the source term vanishes and the equation reduces to the wave equation for (f>, = 0. This can be 

written as an Euler-Poisson-Darboux equation, for which general solutions are available. The full Ernst equation 
has been studied by various methods both in plane symmetry and in the original context of stationary axisymmetric 
space-times; see e.g. the review of Griffiths [3| and references therein. 



VI. LINEARIZED GRAVITATIONAL RADIATION 



To compare with the usual description of linearized gravitational radiation p^], it is convenient to switch tem- 
porarily to Minkowski coordinates {t,x,y,z) defined by \/2x^ — t ± x. Expanding about the Minkowski metric 
r] = diag{— 1, 1,1,1} by 5 = 77-1-/1 consists of expanding about {A, (j>, x, 7) = (1, 0, 0, 0), so one can write A — 1 + a 
and use (a, (/',Xj7) as perturbative fields, each assumed <C 1. Linearizing, the metric perturbation h is given by 

- 2-f{dt^ - dx^) + {a + 2(j))dy'^ + 2{a + 2x)dydz + {a - 2(j))dz^. (44) 

Then the trace of h is 2a + 47 and the trace-reversed metric perturbation h is given by 

a{dt^ ~ dx^) + {2(j) - 2-f)dy^ + 2{a + 2x)dydz + (-20 - 27)^2^. (45) 

Applying the transverse traceless gauge conditions, d°'hab — yields constant a, h^i, — yields a = and ft," = 
yields 7 = 0. Then h = h is indeed transverse: in (y, z) coordinates, 

%)■ 

This verifies the appropriateness of the transverse traceless gauge conditions in plane symmetry. Noting that the 
space-time strain is h/2, this also confirms that (p and x encode the "plus" and "cross" polarizations respectively. 
In the high-frequency approximation, the Isaacson effective energy tensor Q for gravitational waves is defined by 

32TTQab = {dahcddbh-'') (47) 

where the angle brackets denote averaging over several wavelengths [H . Returning to dual- null coordinates, the explicit 
expressions are 

4^e±± = ((9±0)2 + (a±x)'> (48) 

47re+_ = {d+(t)d-<j) + d+xd-x) (49) 
_Le = 0. (50) 

Comparing with ([TTHT5)) . one sees that the radiative components 6±± agree with 0±±, but the other components 
apparently do not. However, this is due to the averaging, as follows. 
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First note that the gravitational wave equation (|4T|l linearizes to the flat-space form 

9+9-$ = (51) 

with general solution 

$ = $+(2;+) + (52) 

as expected. Considering linear superpositions of Fourier modes in the high-frequency approximation, it suffices to 
consider solutions of the form 

$+ = (h+ sin 

+ *X± sm \/2i/±a:=^ (53) 
for constant amplitudes (0±,x±) and angular frequencies {Ltj±,iy±). Then 

uj± cos V2uj±x^ (54) 

d±X = V^Xi'^i cos-\/2i^±a;^ (55) 

and 

47re±± = 20^w^ cos^ V2uj±x^ + 2x\y\ cos^ ^/2v±x^ (56) 

47r±e = 2 cosV2uj+x+ cos\/2a;_x" + x+X-v+v- cos\/2i/+a;+ cos v^x ] 5 (57) 

where 5 = diag{l, 1}. Since (cos^) — \ but (cos) — 0, (-L6) — and similarly 9+_ = 0. Then 

4^(e±±) = 47re±± = 4>lu;l + x\A (58) 

(e+_) = e+_ = (59) 

(±6) = ±e = (60) 

or 

(9) = e (61) 

as expected. Note that the energy densities 6±±/2 have the expected form of squares of amplitudes times angular 
frequencies, with the same numerical factor I/Stt as for electromagnetic radiation in Gaussian units. 

On the other hand, for low-frequency waves, transverse pressure is generally present in _L0 even in the linearized 
approximation, for which O reduces to the energy tensor of a massless complex scalar field in flat space-time: 

STrOab - 29(,$afc)$ - ^ahlf''^c'^^i^. (62) 

The non-zero components (fTT |) -([T3 l) reduce to 

47re±± = (9±(/.)2 + (a±x)' (63) 
47r_Le = (9+09-0 + d+xd-x)^ (64) 

and in particular the transverse shear vanishes, but transverse pressure generally remains. Recall that this is an effect 
for interacting radiation, vanishing for radiation propagating in one direction only. However, if two beams with similar 
amplitude and frequency are passing through one another, the transverse pressure is generally of the same order as 
the energy densities 6±±/2. Although this has been derived here only for plane-symmetric radiation propagating 
in opposite directions, one may expect it to generalize to gravitational radiation from any two sources in different 
directions. 
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